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Abstract 

We use the auxiliary-fields Monte Carlo method for the shell model in the 
complete {pf + 059/2)-shell to calculate level densities. We introduce parity 
projection techniques which enable us to calculate the parity dependence of 
the level density. Results are presented for ^^Fe, where the calculated total 
level density is found to be in remarkable agreement with the experimental 
level density. The parity-projected densities are well described by a back- 
shifted Bethe formula, but with significant dependence of the single-particle 
level-density and backshift parameters on parity. We compare our exact re- 
sults with those of the thermal Hartree-Fock approximation. 
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Nuclear level densities are important for theoretical estimates of nuclear reaction rates 
in nucleosynthesis. The s- and r-pro cesses that involve medium- mass and heavier nuclei 
are determined by the competition between neutron-capture and /5-decay, and the neutron- 
capture cross-sections are strongly affected by the level density around the neutron resonance 
region. Reliable estimates of nuclear abundances often require accurate level densities. For 
example, the abundance of s-process nuclei with non-magic neutron number is (in the local 
approximation) inversely proportional to the neutron-capture cross-section [|l| which in turn 
is proportional to the level density. Most conventional calculations of the nuclear level den- 
sity are based on the Fermi gas model within the grand-canonical ensemble 0. For a gas of 
free nucleons one obtains the well-known Bethe formula. A simple but useful phenomeno- 
logical modification is often adopted, in which the excitation energy is backshifted 0, 
giving a total nuclear level density of 

P^^"(i?.) = g^a-^E^ - A)-te^v/^(^ (1) 

with g = 2. The backshift A originates in pairing correlations and shell effects, while 
the parameter a is determined by the single-particle level-density at the Fermi energy. By 
adjusting the value of a for each nucleus, the backshifted Bethe formula (BBF) (|I]) fits well 
a large volume of experimental data. The value of the parameter, however, is not well 
understood; the Fermi-gas model grossly underestimates the value of a, and cannot account 
for its exact mass and nucleus dependence. Consequently, it is difficult to predict the level 
density to an accuracy better than an order of magnitude. Much less is known about the 
parity-dependence of the level density. The finite-temperature mean-field approximation 
Q offers an improvement over the Fermi gas model but still ignores important two-body 
correlations, especially at low temperatures. 

In this paper we study the nuclear level density in the framework of the interacting shell 
model, in which the two-body correlations are fully taken into account within the model 
space. It should be noted, however, that the finite size of the model space limits the validity 
of such calculations to below a certain excitation energy. The size of the valence shells 
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required to describe the neutron-resonance region for medium-mass and heavier nuclei is 
too large for conventional diagonalization techniques to be practical. However, the recently 
proposed shell model Monte Carlo (SMMC) method makes it possible to calculate thermal 
averages in much larger model spaces by using fluctuating auxiliary-fields. As shown below, 
these methods are particularly suitable for calculations of level densities. 

Nuclei in the iron region play a special role in nucleosynthesis. They are the heaviest 
that can be produced inside normal massive stars, and the starting point of the synthesis of 
heavier nuclei. These nuclei are in the middle of the j9/-shell, and are just beyond the range 
of nuclei where conventional shell model techniques can be applied in a complete p/-shell 
model space Truncated shell calculations [§ were successfully used to describe the 

low-lying states in these nuclei. However, their neutron separation energy, typically ~ 5- 
15 MeV, is too high to justify such truncation. The SMMC method was used to calculate 
thermal properties of ^^Fe in a full p/-shell with the Brown-Richter Hamiltonian. The 
Monte Carlo sign problem of this realistic interaction is overcome through the techniques of 



Ref. |Tn|. However, the statistical errors were too large to obtain accurate level densities. 
Furthermore, the energy range of interest in the iron region {E^ ~ 5-15 MeV) contains 
negative-parity states which are not included in the p/-shell model space. In this letter 
we introduce parity-projection methods for the auxiliary fields, and use the SMMC within 
the full pf- and 05'9/2-shell to calculate total and parity-projected level densities in the iron 
region. This model space is sufficient to describe both positive- and negative-parity states 
for excitation energies up to 20 MeV. To keep the statistical errors small, we construct an 
interaction which is free from the Monte Carlo sign problem, yet realistic enough to describe 
collective features that affect the level density. In particular we present results for ^^Fe, for 
which experimental data are available. 



We adopt an isoscalar Hamiltonian of the form [|TT 

H = Y: ^aUa + 9,P^'^'~^^ ■ P^'''^ - X E ^aO(^-°) ■ 0(^-°) , (2) 

a A 

where 



and (■) denotes a scalar product in both spin and isospin. The modified annihilation operator 
is defined by aj^rn,mt = (— '"^^ ™*«i,-m ,-mt, and a similar definition is used for P^^''^\ To 
conserve the isospin symmetry, the single-particle energies ea are taken to be equal for 
protons and neutrons, and are determined from a Woods-Saxon potential plus spin-orbit 



interaction with the parameters quoted in Ref. |jT2|]. V in (^) is the central part of this 
single-particle potential. The multipole interaction in (^) is obtained (with kx = 1) by 
expanding the separable surface-peaked interaction v{r, r') = —x{dV/dr){dV/dr')6{r — r'). 
The interaction strength x is fixed by a self-consistency relation and we find x = 
0.026 MeV~^fm^ for ^^Fe. In our present calculations we include the quadrupole, octupole 
and hexadecupole terms (A = 2, 3 and 4, respectively). Since our shell-model configuration 
space includes the valence shell alone, core polarization effects are taken into account by 
using renormalization factors kx. We adopt the values /c2 = 2, k^ = 1.5 and k^ = 1, which 
are consistent with a realistic effective interaction in this shell derived by the folded-diagram 
technique |]13|. This interaction satisfies the modified sign rule (suitable for shells with mixed 
parities) |jTO|, and therefore has a Monte Carlo sign of ($) = 1 for even-even nuclei. The 
pairing strength go is determined by using the experimental odd-even mass differences 



for nuclei in the mass region A = 40-80 to estimate the pairing gap. The number-projected 
BCS calculation is then performed for fifteen spherical nuclei with Z = 20, = 28, Z = 28 
OT N = 40 to find the value of qq that will reproduce the estimated pairing gaps. In contrast 
to heavier nuclei |TB|, we find no systematic A-dependence in qq, and a constant mean value 
of go = 0.212 MeV is adopted. 

It is difficult to obtain detailed spectroscopic information on excited states in the 
SMMC method. Instead, it is possible to calculate directly low moments of strength 
functions |jl6[, such as the average energy of the quadrupole excitation i?(Q*^^'°^) = 
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[J:^E^{2t) (2+|g(2,o)|o+) (2+|g(2,o)|o+) ]. Here Q(^'^) is defined by replacing rfV^/rfr 

of O^^'^^ in Eq. by r^. Data on tlie strengtli function of the mass quadrupole mo- 
ment are available from {p,p') experiments in a broad energy range in ^^Fe [|18|. We 
find Eexp(Q(^'°0 ^ 2.16 MeV, in fairly good agreement with the theoretical value of 
^cai(Q^^'°^) = 2.12 ± 0.11 MeV. While this interaction is not fully reahstic (e.g., no spm 
degrees of freedom are included), it seems to reproduce quite well the collective features of 
the nucleus. It is thus reasonable to expect that certain gross properties like the level density 
are described well by this interaction. Note, however, that possible isospin components of 
the effective interaction that push up states of higher isospin are not included in our present 
study since they do not have a good Monte Carlo sign. 

Since both the O/7/2 and OgQ/2 orbits are included in our model space, spurious center- 
of-mass motion may occur. Although this problem has not been fully explored within the 
SMMC framework, it is expected to be unimportant in our case. The energy difference 
between O/7/2 and OgQ/2 is 9.6 MeV, so we expect spurious states to appear at excitation 
energies around 10 MeV and higher. However, their density is comparable to that of the 
non-spurious states but at about 10 MeV lower in energy, and is thus a negligible fraction 
of the total density at the actual excitation energy. 

In the SMMC the energy is calculated as a function of inverse temperature f3, from the 
canonical expectation value of the Hamiltonian E{f3) = {H)p through an exact particle- 



number projection |19| of both protons and neutrons. The partition function Z[j3) is then 
determined by a numerical integration of E{(3) 

ln[Z{f3)/Z{0)] = - f'dp'EiP'), (4) 
Jo 

where Z{0) = Tr 1 is just the total number of states within the model space. The level 
density p{E) is the inverse Laplace transform of Z{P), and is calculated in the saddle-point 
approximation from 

p{E) = (27r/3-2^)-i/2g5 . 

S{E) =(3E + In Z{(3) , (3-^C{f3) = -dE/dp . (5) 



Here j3 = (3{E) is determined by inverting the relation E = E{f3), and C is the heat capacity 
calculated by numerical differentiation of E{f3). 

We now introduce parity-projection techniques in the SMMC to calculate parity- 



projected observables (see Ref. [|17| for parity-projected thermal mean-field approxima- 
tion). Using the Hubbard- Stratonovich representation for e~^^ and the projection operators 
P± = (1 ± P)/2 (P is the parity operator) on states with positive and negative parity, re- 
spectively, we can write the projected energies E±[(3) = Ti[HP±e~^^) /Ti{P±e~^^) in the 
form 

^^^^ SD[a]W{a)[l±Ua)/a<y)] ' 

The integration over the auxiliary fields a is performed with the usual Monte Carlo weight 
function W{a) = G(o")C(cr), where G is a gaussian factor and ({a) = Tr Ua is the par- 
tition function of the non-interacting propagator Ua p. In (H) Cp{cr) = Tr(Pf/o-) and 
{H)pa- = Tt{HPU„)/Tt{PU„). Since the parity operator can be expressed as a product of 
the corresponding parity operators of each of the particles, it follows that the operator Pf/^- 
can be represented in the single-particle space by the matrix VU^j, where is the matrix 
representing Ua in the single-particle space, and P is a diagonal matrix with elements (— )^' 
{ii is the orbital angular momentum of the single-particle orbit i). This representation al- 
lows the calculation of Cp(^) ^"^^ {H)pa through matrix algebra in the single-particle space, 
similar to the calculation of ({a) and {H)^., except that the matrix U„ is replaced by VU„. 
Once we calculate E±{P), we can proceed to calculate the densities p±{E) as for the total 
level density. 

In the following we present results for ^^Fe. For each f3 we used a Monte Carlo time 
slice of A/5 = 0.03125 MeV~^, and collected more than 4,000 samples. To describe the level 
density as a function of the excitation energy E^, we also need to know the ground-state 
energy. The latter is calculated by extrapolating E{P) to /3 — cxd. Fig. |1] shows the energy 
E as a function of (3. The SMMC results are the solid squares and include statistical errors 
(although the errors are often too small to be visible in the figure). We compare our results 
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with those of the thermal Hartree-Fock approximation (HFA) where we observe large 
deviations at low temperatures. In the HFA, a shape phase-transition occurs around /5 ~ 1.3 
MeV~^ from a spherical configuration (at higher temperatures) to a deformed one, and its 
signature is observed in the bending of E{(3). The inset to Fig. |I| presents E± as a function 
of /3, calculated using the parity-projection technique of Eq. (j^). Because of the energy gap 
between pf and 0(^9/2, E^{P) is notably higher than E^{[3) at low temperatures. 

The SMMC and HFA entropies S (versus E) and heat capacities C (versus (3) are shown 
in Fig. ^. We observe a substantial enhancement in the SMMC entropy over the HFA 
entropy, which is due to the full inclusion of the two-body correlations. The heat capacity 
is useful for determining the range of excitation energies for which the present model space 
is sufficient. The decrease of C at high excitation energy (i.e., small (3) is associated with 
the truncation of the single-particle space, and we conclude that our present calculation is 
meaningful up to Ex ~ 20 MeV. The discontinuity of the heat capacity at /? ~ 1.3 MeV~^ in 
the HFA is a signature of the shape transition, but this effect is washed out in the SMMC. 

The SMMC total level density is shown in Fig. ^ (left panel) as a function of E^- Although 
it is difficult to measure the total level density directly, it can be reconstructed from a 
few parameters that are determined experimentally. The solid line in Fig. ^ shows this 
experimental level density with the BBF parameters of a = 5.80 MeV~^ and A = 1.38 MeV 
||20|| . Our SMMC result is in excellent agreement with the experimental level density. The 
slight discrepancy at low energies may be ascribed to deviations of the moment-of-inertia 
parameter from its rigid-body value: a rigid-body moment was assumed in deriving the 



the experimental values of a and A. We can also use our microscopically calculated level 
densities to extract the level density parameters via a fit to Eq. (|lD . Using the energy range 
4.5 MeV < Ex < 20 MeV, we obtain a = 5.780 ±0.055 MeV"^ and A = 1.560 ±0.161 MeV. 
We note that the statistical errors of the present calculations are substantially smaller than 
those of Ref. [Q, where thermal properties of ^^Fe were calculated using a realistic interaction 
in a smaller configuration space (p /-shell). Consequently, accurate level densities can be 
calculated in the present work. Also shown in Fig. |^ is the level density in the HFA, 



where the excitation energy has been corrected by the difference between the mean-field 
and SMMC ground state energies. The SMMC level density is significantly enhanced in 
comparison with the HFA level density. The kink around 9 MeV in the HFA level density is 
related to the shape transition, and disappears in the SMMC. Although the Hartree-Fock- 
Bogoliubov approximation is expected to improve the HFA because of the T = 1 pairing 
component, the shortcomings of the mean- field approximation are likely to remain. Indeed, 
even when the pairing force is omitted from the interaction, we still find that the HFA level 
density has a kink and deviates strongly from the exact SMMC level density. 

The Fermi gas model predicts equal positive- and negative-parity level densities at all 
energies. However, this is unrealistic in the neutron resonance regime where the neutron 
resonance energy is comparable to or even smaller than the energy gap among major shells. 
The SMMC results for the parity-projected level densities of ^^Fe are shown in the right 
panel of Fig. |[ They can be well-fitted to a BBF (|l|) with g = 1, but with parity-specific 
parameters a± and A±. We find a+ = 5.611 ± 0.073 MeV"\ A+ = 0.550 ± 0.196 MeV 
and a_ = 6.209 ± 0.625 MeV"\ A_ = 3.172 ± 1.637 MeV. We remark that negative-parity 
states in ^^Fe are possible only when the 5^9/2 level is populated. Because of the energy gap 
between the pf and gg/2 orbits we expect the negative-parity level density to be lower than 
the positive-parity density at low energies. Thus the backshift A_ should be larger than A+, 
in agreement with our results. Both Aj. are significantly different from A of the total level 
density. On the other hand, a_|_ is rather close to a, while a_ is larger than a_|_ (and a). At 
high excitation energies the Fermi-gas model is expected to be a reasonable approximation, 
implying the approximate equality of positive- and negative-parity level densities (in our 
case p+ ~ p_ above ~ 17 MeV). Therefore, in the low energy region the negative-parity 
density is expected to rise more quickly as a function of energy, i.e., a_ > a_(_. This relation 
is confirmed by the present calculations. So far there has been no systematic study of the 
parity-dependence of the level density parameters. It would be interesting to investigate 
how this parity-dependence affects the neutron-capture reaction rates. 

In conclusion, we have used the auxiliary-field Monte Carlo methods to calculate the 
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level density of ^^Fe in the complete pf- and g'9/2-sliell, and found remarkable agreement 
with the experimental level density. The SMMC calculations are an important improvement 
over the finite temperature Hartree-Fock approximation. We have introduced a novel parity- 
projection technique in the SMMC, which allows us to study the parity dependence of both 
the single-particle level density parameter a and the backshift parameter. Work in progress 
includes a systematic study of level densities for nuclei in the {pf + gi9/2)-shell. 
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FIG. 1. The total energy £^ as a function of P for ^°Fe. The SMMC values are shown by solid 
squares, while the HFA values by a dashed line. The inset shows the positive (circles) and negative 
(triangles) parity-projected SMMC energies. 
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FIG. 2. Right: the entropy as a function of energy E. Left: the heat capacity C as a function 
of p. The conventions are as in Fig. 1. 
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FIG. 3. Level densities of ^"^Fe. Left: total level density. The SMMC level density (solid 
squares) is compared with the HFA level density (dashed line). The solid line is the experimental 



level density |20|. Right: positive- and negative-parity level densities in the SMMC. The conven- 
tions are as in the inset to Fig. 1. The dotted lines are the fit to Eq. (||) with the parameters 
quoted in the text. 
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